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Abstract. In this paper, multi-component generalizations to the Camassa-Holm equation, 
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1 Introduction 

Integrable systems solved by the inverse scattering method usually arise from shallow water 
wave, physics, optical communication and applied sciences etc. Integrable systems have many 
interesting properties, such as Lax-pair, infinite number of conservation laws and Lie-Backlund 
symmetries, multi-solitons, Backlund transformations and bi-Hamiltonian structure etc. [1,44], 
which are helpful to explore other properties of integrable systems [1,44,56]. 

It is of great interest to study geometric aspects of integrable systems. So far, very few inte- 
grable systems were found to have geometric formulations. The relationship between completely 
integrable systems and the finite-dimensional differential geometry of curves has been studied 
extensively. It turns out that some integrable systems arise from invariant curve flows in certain 
geometries [2-7,9-15,21,23,24,26,28-42,45,47-52,55,57,59,60]. The pioneering work on this 
topic was done by Hasimoto [24] . He showed that the integrable nonlinear Schrodinger equation 
(NLS) 

i(t)t + (t)ss + \(t)?(t> = ^ 
is equivalent to the system for the curvature k and r of curves 7 in 

^sss ^ s ^ ss / \ 

Kt = -2tKs - KTs, Tt = 5 2tTs + KKs (1) 



*This paper is a contribution to the Special Issue "Symmetries of Differential Equations: Frames, Invariants 
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via the so-called Hasimoto transformation (p = Kexp{i T{t, z)dz). Indeed, the system (1) is 
equivalent to the vortex filament equation 

7t = 7s X Iss = Kh, (2) 

where b is the binormal vector of 7. Mari Beffa, Sanders and Wang [39, 51] noticed that 
Hasimoto transformation is a gauge transformation relating the Frenet frame {t,n, b} to the 
parallel frame {t^,n^,b^}. It is also a Poisson map which takes Hamiltonian structure of the 
NLS equation to that of the vertex filament flow [29]. The Hasimoto transformation has been 
generalized in [51] to the Riemannian manifold with constant curvature, which is used to obtain 
the corresponding integrable equations associated with the invariant non-stretching curve flows. 
The parallel frames and other kinds of frames are also used to derive bi-Hamiltonian operators 
and associated hierarchies of multi-component soliton equations from non-stretching curve flows 
on Lie group manifolds [3,4,31,39]. The KdV equation, the modified KdV equation, the Sawada- 
Kotera equation and the Kaup-Kuperschmidt equation were shown to arise from the invariant 
curve flows respectively in centro-equiaffine geometry [7,9,48], Euclidean geometry [21], special 
affine geometry [11,35] and projective geometries [11,30,41]. 

The integrable systems with non-smooth solitary waves have drawn much attention in the 
last two decades because of their remarkable properties. The celebrated Camassa-Holm (CH) 
equation 

mt + urrix + 2uxm + aUx = 0, m = u — Uxx, (3) 

was proposed as a model for the unidirectional propagation of the shallow water waves over 
a flat bottom, with u{x,t) representing the water's free surface in non-dimensional variables [8]. 
It was also found using the method of recursion operators by Fokas and Fuchssteiner [19] as 
a bi-Hamiltonian equation with an infinite number of conserved functionals. Geometrically, 
the Camassa-Holm equation arises from a non-stretching invariant planar curve flow in the 
centro-equiaffine geometry [9], and the periodic CH equation (3) describes geodesic flows on 
diff (S^ X M) with respect to right-invariant Sobolev metric for a = [13, 14, 28] and Bott- 
Virasoro algebra for a 7^ [40] . A dual version of the Ito system is the two-component Camassa- 
Holm equation [46], the periodic two-component CH equation also describe geodesic flows on 
an extended Bott-Virasoro algebra [22]. 

It is remarked that all nonlinear terms in the CH equation are quadratic. In contrast to 
the integrable modified KdV equation with a cubic nonlinearity, it is of great interest to find 
integrable CH-type equations with cubic or higher-order nonlinearity and non-smooth solitary 
waves. To the best of our knowledge, two scalar integrable CH-type equations with cubic 
nonlinearity have been discovered. The first equation is [18,46,49] 

nT-t + ^ {{u^ + Sul)m) ^ = 0, m = u + 6uxx, (4) 

where 5 = ±1, and the second one is the so-called Novikov equation [25,43] 

mt + Suuxin + u nix = 0, m = u — Uxx, 

which are completely integrable, and admit peaked solitons. Recently, systems of CH-type 
equations with cubic nonlinearity were also obained [20,55]. 

The CH equation can also be derived by the tri-Hamiltonian duality approach basing on 
bi-Hamiltonian structure of the KdV equation. Other examples of dual integrable systems 
obtained using the method of tri-Hamiltonian duality can be found in [18,46]. Nonlinear dual 
integrable systems, such as the CH equation and the modified CH equations, are endowed with 
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nonlinear dispersion, which in most cases, enables these systems to support non-smooth soliton- 
like structures. It was remarked in [23] that the modified CH equation (4) can be regarded 
as a Euclidean-invariant version of the CH equation (3), just as the modified KdV equation is 
a Euclidean-invariant counterpart to the KdV equation from the viewpoint of curve flows in 
Klein geometries [9,10,21,48]. 

The aim of this paper is to provide geometric formulations to multi-component integrable 
systems admitting non-smooth solitons. We shall show that several multi-component integrable 
systems with non-smooth solitons, such as a dual version of the Schrodinger equation [17,46], 
the complex CH equation and multi-component modified CH equations arise from the invariant 
curve flows respectively in three-dimensional Euclidean geometry, Mobius sphere and the n- 
dimensional unit sphere S"(l). To obtain integrable systems relating to these geometric flows, 
we shall use the scale limit technique. The outline of this paper is as follows. In Section 2, 
a non-stretching invariant binormal curve flow in is introduced and studied. Making use of 
the system for curvature and torsion corresponding to this flow, we obtain a novel integrable 
Schrodinger equation by a scale limit approach, which is completely integrable system and can be 
obtained by the so-called tri-Hamiltonian duality approach [18,46]. In Section 3, we give a brief 
discussion on Mobius 2-sphere P0(3, 1)/H and the n-dimensional sphere SO(n-|- 1)/S0(n), the 
Cartan structure equations for curves in both geometries are reviewed, which will be used in 
subsequent sections to study curve flows in both geometries. In Section 4, we consider the 
non-stretching curve flows in Mobius 2-sphere. It is shown that the complex Camassa-Holm 
equation and complex Hunter-Saxton equation describe the non-stretching curve flows in Mobius 
2-sphere. The bi-Hamiltonian structure for the complex Camassa-Holm equation is obtained. In 
Section 5, we study non-stretching curve flows in the n-dimensional sphere S"'(l). Interestingly, 
we find that a multi-component modified CH equation (a multi-component generalization of 
the modified Camassa-Holm equation) is equivalent to a non-stretching curve flow in §"'(1). 
Integrability of the system is identified. Finally, Section 6 contains concluding remarks on this 
work. 



2 An integrable nonlinear Schrodinger equation 

We consider the flows of space curves in M'^, given by 

-ft = Un + Vh + Wt, (5) 

where t, n and b denote the tangent, normal and binormal vectors of the curves, respectively. 
The velocities U, V and W depend on the curvature and torsion as well as their derivatives with 
respect to arc-length parameter s. The arc-length parameter s is defined implicitly by ds = hdp, 
h = \'y'{p)\, where p is a free parameter and is independent of time. We denote by k and r the 
curvature and torsion of the curves, respectively. Governed by the flow (5), time evolutions of 
those geometric invariants are given by [24,42] 

^ - tV + kw] n+ +tU ]h, 



'dU 

iv=-[—-TV + KW\t + 
OS 



1 d (dv ^\ T (du 

KOS \ OS J K \ OS 



b, 



1 d fdV , \ r fdU 

{ ^ + tU] + - { tV + kW 



K.ds yds J K, \ds 



n. 



h = 2h{^-.U] (6) 



4 



C.Z. Qu, J.F. Song and R.X. Yao 



and 



dr d 
dt ds 
dK d'^U 



1 d (dv ^\ T (du ^\ r , 

KOS \ OS / K \ OS 



KtU + K 



dV 
ds ' 



dt ds"^ 



ds 



\uds'-2r^-^V. 

OS OS 



(7) 



Assuming that the flow is intrinsic, namely the arc-length does not depend on time, it implies 
from (6) that 



Ws = kU. 

In terms of (7), one finds that the complex function 

i / T{s',t)ds' 



(8) 



K7], 



r] = exp 



satisfies the equation [24,42] 



[ds^ 



ds 



ds'cj) ] (Ut]) 



ds^ 



+ [ i— + + (I) I ds'r^-i^ I ds'^]{VT]), 



(9) 



where (p denotes the complex conjugate of (j). 

Let U = 0, V = K, where k is a real function, then (8) implies that W = Ci, where Ci is 
a constant. Setting Ci = 0, we derive from (9) the celebrated Schrodinger equation 



i(l)t + (t>ss + ^\^\'^(t> = 0. 



(10) 



Let U = —Kg, V = —KT, then W = + C2, where C2 is a constant. Letting C2 = 0, we 

find that (j) satisfies the mKdV system 

3 

(l)t + (t)sss + -\(k\^(l)s = Q. 

In the following, we shall consider the case U = W = Denote 0{s,t) = r(s', t)ds', 
g = Vrj. It follows from (9) that (p satisfies the equation 



+ dss + \(p\'^g — 4> J g{cos 9 — i sin 6)kg/ds' = 0. 



(11) 



Set u = KcosO, V = Ksm9, g = gi+ig2, then the equation (11) is separated to two equations 

ut = -g2,ss -vd~'^[K{gi cos 6* + 52 sin6')s], 
vt = gi,ss + ud~^[n{gi cos 6* + 52 sin6')s]. 

Furthermore, letting u = u + Vs, v = v — Us, and choosing the binormal velocity V to be 
V = dj^l^u^ + v'^^ ^/ k], we find that u and v satisfy the system 



(n + Vs)t = -g2,ss -{v - Us) {u^ + v^) , 

(v - Us)t = gi,ss + {u^ + v"^) {u + Vs) 



(12) 
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{U + Vs^ + {V - UsY ' [U + VsY + {V - UsY ' 

Applying the scaling transformations 

s I — > s, 1 1 — > e^t, u I — > e^^u, v i — > e~^v 
to system (12) produces 

(n + Vs)t = -€^g2,ss - {v- Us) {v? + v^) , 

{v - Us)t = e^9i,ss + {u^ + w^) (u + Vs). (13) 
Expanding u and v in powers of the small parameter e 

u{t, s) = uo{t, s) + eui{t, s) + e^U2{t, s) + ■ ■ ■ , 
v{t, s) = vo{t, s) + evi{t, s) + e^V2{t, s) + • • • , 

and plugging them into system (13), we find that the leading order terms no(f,s) and fo(t,s) 
satisfy the system 

(no + vo^s)t + {ul + vl) {vq - Uq^s) = 0, 

(^0 - uo^s)t - {ul + vl) (uo + Vq^s) = 0. (14) 

Again we use the notation (t){t, s) = uo{t, s) + ivo{t, s), then it is inferred from (14) that </>(t, s) 
satisfies the equation 

i{(Pt-icj)ts) + \cl)W(l)-i<Ps)=0, (15) 

which is a dual version of the Schrodinger equation (10), and can be obtained by the approach of 
tri-Hamiltonian duality [18,46]. Equation (15) is formally completely integrable since it admits 
bi-Hamiltonian structure [46] 

op op 
where p = (p — i(f)s, Si and £2 defined by 

£i=ds + i and £2 (F) = pd~^ {pF - pF) 
are compatible Hamiltonian operators, while 

Hi= I (|0|'-i#s)ds= / p^ds, H2 = \ I |0|2(|</.|2-i#,)ds. 

are the corresponding Hamiltonian functionals. 



3 Mobius sphere PO(3, 1) /H and unit sphere SO(n + 1) /SO(n) 

In this section, we give a brief account of Mobius 2-sphere P0(3, 1)/H and unit sphere ^"(l) = 
SO(n + 1)/S0(n). Please refer to the book [53] for the details of the two geometries. 
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3.1 Mobius 2-sphere 

Let {uo,ui,U2,us) G M'^, we define the inner product on by 

{x,y) = x^As^iy, 
where x,y E M^, and the matrix A3 1 is 



A3,l 



/ -1\ 

10 

10 

V-1 / 



A vector field x G M"^ is said to be light-like, if it satisfies {x, x) = 0. All the light-like vector 
fields form a set L, which is called optical cone, defined by the equation 

2x0x3 — Xi — X2 = 0. 

Clearly it is homogeneous, namely for any A G M, if x G L, then Ax G L. The projectivisation 
of L is said to be Mobius 2-sphere, which is isomorphic to . Recall that 

0(3,1) = {ge GL(4,R) : /A3,i5 = Agj} , 

and the Mobius group is defined to be P0(3, 1) = 0(3, 1)/±I. We denote 

H = {h£FO{S, 1) : /iN = [63]} 
0\ 

A G 0(3,1); AG 0(2), a GM+, 'yG™^ 
f]'^ a) 

where [63] denotes the equivalent class of 63 in P(M*^), [63] = (0,0,0,*). A straightforward 
computation gives 

am a 
-V V = -\ 
2 2' 




-iT 7 CI *^ II ii2 

r] = aA V, = -V v = -\\v\\ . 



It is easy to verify that the group P0(3, 1) acts on the Mobius sphere transitively (the group 
action is the usual conformal transformation). 

Definition 1. The Klein pair (P0(3, 1),H) is called the Mobius 2-sphere. 

For any g G P0(3, 1), there exists a unique decomposition around the identity of the group 
1 

g = 51505-1 

where a G M+, A G 0(2), u,v £ R^, h = gigo G H, u in g^i part represent a local coordinate of 
the point in P0(3, 1)/H. For the corresponding Lie algebra g, there exists also a decomposition 

0-1, 





-€ 0\ 

5 Gflo, I q\ Gg_i, 
e/ 

p,q G M?, e G M, S G 0(2). The Lie algebra of the isotropy group H is I) = Qi (B go while 
g/f) = 0_i is identified to the tangent space of the conformal sphere P0(3, 1)/H. 
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3.2 n-dimensional sphere §"(1) = SO(n + l)/SO(n) 

The n-dimensional unit-sphere is also a homogeneous space M = G/H = SO(n + 1)/S0(n). 
The corresponding Lie algebra has the following Cartan-Killing decomposition 

so(n + 1) = f) e m = so(ra) © M'', 

with 

where p G M", G so(n), and the decomposition satisfies 

[f),f)]cf), [{),m]cm, [m,m]c{), 

where m is identified to the tangent space T^M = M" of M = S0(n + 1)/S0(n). The flat Cartan 
connection of principle S0(?7.) bundle SO(n + 1) §" is given by the Maurer-Cartan form of 
Lie group SO(n -|- 1). The Cartan structure equation reads as 

Q = duj + -[uj,uj] = 0, 
where 50 (n -|- l)-valued one-form lo is decomposed to 

where g/f)-valued 9 represents a linear coframe, f)-valued ujh represents a linear connection on S". 
The corresponding Cartan structure equation is separated to 

^ = de + ^[uH,9] + ^[e,ujH] = o 

and 

where Z and 9K are called torsion and curvature forms, respectively. 

4 Curve flows in PO(3, 1)/H and the complex CH equation 

For the Mobius geometry P0(3, !)/-&, its Cartan connection takes values on g = o(3, 1), with 
the form 

\ T]^ ej 
and the Cartan structure equation reads as 

Q = doj + -[uj , Lo] = 0. 

Let's consider the invariant curve flows for curves 7(x, t) = {ti.i{x, t), %L2{x, t)) on the conformal 
sphere = P0(3, l)/H, where x denotes the parameter of the curves, t is the time variable, 
ui and U2 denote the local coordinates on . Let 7t = 7* ^ be the evolutionary vector field of 
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the curves, — 1*'§~ denotes the tangent vector of the curves. Assuming that the curve flow 
is intrinsic, namely the parameter x for the curve does not depend on time t, we have 

bx,it\ = 0. 

It was shown in [16, 36] that there exists a conformally equivariant moving frame (the Frenet 
frame) p = p{x,t) S P0(3, 1) along the curve 7(x,t) C M^. Let Dx and Df denote respectively 
the vector field ^ and ^ along the curves p in P0(3,l), then the Frenet formulae for the 
conformally parametric curves is 

Px = pu}{Dx), 

with 



(^{Dx 





1 







h 








1 













Vo 


ki 


k2 


0/ 



(16) 



where ki and k2 are the conformally differential invariants for the curve 7(3;, i). The time 
evolution for the frame p{x, t) can be written as 



pt = pio{Dt 



where 



w(A 





hi 


h2 




fi 





—a 


hi 


f2 


a 





h2 


VO 


h 


f2 





(17) 



and e, a, fi, hi (i = 1,2) are some conformal differential invariants related to ki and k2, to be 
determined. By the Cartan structure equation, one gets 



n{Dx,Dt 



dt 



CbiDx) 



[CbiDx),u;{Dt)] = 0. 



Plugging (16) and (17) into (18) results in the following equations 

e = -hi^x, OL = h2,x, 
fi= ex + kihi + k2h2, 
f2 = Q-x + k2hi - kih2, 
ki,t = fi,x - kie + ak2, 
k2,t = f2,x - ek2 - aki. 



(18) 



(19) 
(20) 
(21) 
(22) 
(23) 



where (19) is the torsion part (i.e., the g_i part) of the Cartan structure equation (18), which 
can be written as 



-dx 
dx 



= Ji 



Inserting (19) into (20) and (21) gives 

fi = -hi^xx + kihi + /C2/12, /2 = h2,xx + k2hi - kih2. 



(24) 
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Substituting (19) and (24) into (22) and (23), we obtain the evolution equations for the curva- 
tures ki and k2 [32,36] 

ki,t = -hi,xxx + 2A;i/ii,a; + ki^xhi + 2k2h2,x + k2,xh2, 
k2,t = h2,xxx - 2A;i/i2,x - ki^xh2 + 2A;2/ii,x + ^2,z^i, 
which is equivalent to 

ki\ _ f-d^ + kid + dki k2d + dk2 \ (hi 
k2)~\ k2d + dk2 -kid-dkj [h2 



= J 2 



n 



The following cases are considered. 

Case 1. Setting ki = 1/2 + m = 1/2 + n — Uxx, ^2 = 
in (25), we obtain the new two-component CH equation 

mt + 2ux'm + urux — 2vxn — vux = 0, 
nt + 2uxn + uux + 2vxTn + vrux = 0. 

The above system admits the following Lax-pair 

4>x = U(j), 4>t = V(t>, 

with 

/O I + Xm Xn \ 

1 l + Xm 
An 
\0 1 / 

where fi = ^(A"^ + n) + X{vn — um), f2 

Case 2. Setting ki = m = —Uxx, k2 = n = —Vxx, hi 
complex Hunter-Saxton equation 

mt + 2uxm + urUx — 2vxn — vUx = 0, 
nt + 2uxn + uUx + 2vxTn + vrrix = 0, 

which admits the Lax-pair (27) with 



U 



(25) 

Vxx, hi = 1 — u and h2 = v 

(26) 
(27) 



U 



and 



/ 


-Ux 


/l 


/2 


°\ 


A- 









/l 




V 







/2 


v 





A-1 - w 







— X{un + vm). 



1 — u, h2 = V, we arrive at the 



and 





AlixX 


-Au^.^ 


\ 








1 








'XUxx 

















->^Vxx 








Vo 


1 













/ 


-Ux 


X{vn — um] 


—X{un + f m) 





\ 


A^ 











A(vn — um) 






V 


-Vx 







—X{un + 




v 





X-^-u 




V 







V 



Case 3. In (25), letting hi and /12 satisfy 

hih2 



ki 



1 hi -hi 



'^{hl + hlf ^ {hl + hiY 

we get the two-component Harry-Dym equation 

hi^t = {h\ — Shih^hi^xxx + {hi — 3hlh2)h2,xxx, 
h2,t = {h2 — 3/1^/12) /il.xra — (^1 — 3hihl)h2,xxx- 



k2 



v2 ' 
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Remark 1. Mari Beffa [32,36] showed that the complex KdV equation arises from the invariant 
curve motion in Mobius 2-sphere. Indeed, taking hi = —ki, h2 = —k2 in (25) yields the complex 
KdV equation 



2,t 



'^2,xxx 



3ki k 



lK2,x 



3k2k 



Its Hamiltonian structure J2, see (25), was originally derived in [32]. One can see that the bi- 
Hamiltonian structure Ji and J2 of the complex KdV equation comes from the Cartan structure 
equation for the conformal invariant curve flow. According to the decomposition of the Lie 
algebra 

g = 01 ©So ©0-1, 

the Cartan curvature form Q is decomposed to = l^i + Oo + ^-i; where Ji comes from torsion 
part of the structure equation, i.e., the 0_i part, and J2 arises from the r^i part. In the sequel, 
we will show that the complex CH equation admits a bi-Hamiltonian structure J7i and ^2- It 
turns out that the complex CH equation is a dual version of the complex KdV equation (in the 
sense of [46]). 

It is well-known that the CH equation is a bi-Hamiltonian system [8] 



mt = J 



6H1 
6m 



V 



SH2 
5m 



where = —{md + dm) and T> = —{d — d^) are Hamiltonian operators. Hi = — 2 / ^^dx and 
H2 = ~5 / u^u^+u^^dx are the corresponding Hamiltonian functionals. As for the CH equation, 
the complex CH equation can be obtained by the approach of tri-Hamiltonian duality [46]. 
Indeed, we have the following result. 

Theorem 1. The complex CH equation (26) is a bi-Hamiltonian system, which can be written 
as 



Ji 



(5H2\ f6Hi\ 

= J2 



5m 
6H2 
\ 6n J 



6m 
6H1 



with 



and 



Ji 



Hi 



93 _ a 

d-d^ 



[vn — um)dx, 



J2 



md + dm nd + dn 
nd + dn —md — dm 



= 2 y [''^('"^ + ^z) ~ m(3w^ -|- w^) — 2vUxVx] dx. 



Proof. Clearly, J'l and J'2 are skew symmetric. To prove they are Hamiltonian operators, it 
suffices to prove that the Poisson bracket defined by J2 satisfies the Jacobi identity. 
The bi-vector associated with J2 is defined by [44] 

1 



e 



J2 



(t9 a J2'd)dx 



{mO A Ox + n9 A Cx + nC A 9x - mC A (x) , 



where = (^, C); ^ C denote the basic unit vectors corresponding to m and n, respectively, 
the notation A denotes the usual inner product between ■& and J^2'&- It suffices to show that the 
Schouten bracket vanishes, namely \J2-,J2\ = 0. In terms of 



pr V 



m 



2m0x + mxO + 2nCx + rixC, P^^ ^iai?^") ^ ^""^^ + ""^^ ~ ^"^^^ ~ 
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a direct computation shows 

(2nCx A 6* A 6*^ + n^C A 6* A 6*^ + 2ne^ AOaC^- m^C A 6^ A Cx 

+ n^r^' A C A Or, - 2mCx A C A 6*^ - 2mBx ^ C ^ Cx - rrixO A C A ex)dx = 0, 

where the skew-symmetric property for the wedge product is used. 

Next, we prove that the Hamiltonian operators Ji and ^2 are compatible, i.e. 

[Ji , J2] + [J2Ji] = wvj^^{ej.) + v^v j,^ {QjJ = 0. 

Note that 

piVj^^iQ jj = 0, picvj^^{m) = Oxxx - Ox, W^j-^i}{n) = Cx - Cxxx- 
Through integration by parts, we get 

[J1J2] + [J2,Ji] =prvj^^{ejJ+pTVj^^{ejJ 

^xxx /\0 A9x — Cxxx /\ A Cx + Cx ^ C ^ Ox — Cxxx A C A ^a; 
:xx 

A C A Cx + A C A Cx)dx 
( Cxxx A 6 A Cx Cxxx A C ^ Ox Oxxx A C A Cx)dx 

{-Cxxx AOACx- Cxxx /\ C ^ Ox - Ox A Cx ^ Cxx - Ox A C A Cxxx)dx 
{-Cxxx AOACx-OxACx/\ Cxx)dx = 0. 

Thus J'l and ^2 are a Hamiltonian-pair. Let's write the complex CH equation (26) as follows 
/6Hi\ 



J2 



5m 
6H1 



It is easy to get 



5Hi 6H1 



-u, 



dm 6n 
Thus 

6H1 5Hi 

— — = — m, -z — = n. 
ou ov 

Hence we deduce that 

1 f 

Hi = 2 / ~ '^™^dx. 
To compute fl2, we write (26) as 



mt = -d\ ^li^ - \^u\ - uuxx - ^v'^ + ^vl + vvxx / , 
nt = d {uvxx + vuxx + UxVx — 'iuv) . 
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Since H2 satisfies 

mt = -9(1 - 9 — = -o^—, nt = 5(1 - o^)^— = d^—. 

dm ou on ov 

It follows that 

SH2 

— = UVxx + VUxx + UxVx — OUV. 

dv 

Note that 

^2 ^2 ^ ^ f / 2\ ^ 

-U - - UUxx = / [u +uux)dx, 

6_ 
6u 
6_ 
6v 
6_ 



2 J 2 



uv^dx = 2uv, 



Uvl - UxVVx )dx= -v1+ VVx 



Hence we arrive at 

H2 = - I [u + — ?)UV — uVx — 2uxVVx) dx. ■ 

5 Curve flows on S"(l) and multi-component 
modified CH equations 

Assume that 7(x,t) is a curve flow on unit sphere §"'(1) = SO(n + 1)/S0(n), which satisfies 
II7II = 1, where x is the invariant arc- length parameter, t is the time. The natural frame of 
the curve 7 S S"(l) is {ei = 7^;, 62, ... , e„}. Let /? = (eo = 7, ei, . . . , e„) G SO(n + 1) be the 
lift from S"(l) to bundle space SO(n + 1), and Dx and Dt denote respectively the tangent and 
evolutionary vector field. It follows that 

Px = pCj{Dx), (28) 
where uj is the Cartan connection 



'^{Dx 



/O -1 

1 -fc'^ I , O G soin - 1) 



k 



(the natural frame formulae (28) for curves on the sphere comes out from the Frenet formulae [15] 
through the Hasimoto transformation). Here k = {ki,k2, ■ ■ ■ ,kn-i) is the natural curvature 
vector of 7. 

Assume that the curve flow is governed by 

It = fei + hie2 + /i2e3 H h K-ien, 



where the tangent velocity / and normal velocities hi (z = 1, 2, . . . , n — 1) depend on the curva- 
tures and their derivatives with respect to arc-length x. 
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The induced time evolution for the frame is 

pt = puj{Dt), 

with 



f 

\h e 



9 G so(n - 1), 

e & J 

^ is a unknown vector, which will be determined later by the structure 



where h,^ G 
equations. 

First, we assume that the flow is intrinsic, namely, the distribution {Dx,Dt} satisfies 
[Dx,Dt] = so that the integral submanifold is a smooth two-dimensional surface on Lie group 
SO(n + 1). Making use of the Cartan structure equation 



one gets the following equations 

fx = {k, h), 

i= hx + fk, 

Qx = k<g)i-i(g) k, 

kt = - Qk + h, 



0, 



(29) 
(30) 
(31) 
(32) 



where (29) is the arc-length preserving condition. 

For convenience, the following notations are used. For any a, b £ W'-^, {a,b) denotes the 
usual Euclidean inner product, i.e., (a, 6) = a^b, a0 5 denotes the tensor product, namely 



/ aibi 



aib2 



aiK-i \ 
a2bn-i 

an-iK-iJ 



(33) 
(34) 



Define aAb = a(^b — b0a. From (29) and (30), it follows that 

/ = d-\k, h), i= hx + id-\k, h))k. 

In view of (31), we have 
Q = d-HkAi). 

Plugging (34) and (33) into (32) leads to the equation for the curvature vector 

h = hxx + {k, k)h+ {dx\k, h))kx + (9-1(4 A h))k + h, (35) 

where the identity for vectors 

{a Ab) • c = {b,c)a — {a,c)b 

was used. 

Analogous to the derivation for the CH equation [18] and the modified CH equation [46], we 
restrict our attention to the following cases. 
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Case 1. h = Ux, k = rh = u + Uxx- In this case, the tangent velocity / is determined by 

/ = d^^{u + Uxx,Ux) = ^(I'wp + l^^xP) +co = ^Q + co, (36) 

where cq is an integration constant. Substituting (36) with cq = —1 together with the expressions 
for k and h into (35), and noting that 

{d~^{kx Ah))k= [d^^{uxxx A Ux)) m = (uxx A Ux) m 

= — {u A Ux) rh — {ux Am)rh = — {u A Ux) m — {ifi, rh)ux + -Qx^t-, 

we arrive at the multi-component modified CH equation 
1 

= - {rnQ)^ - {u A Ux)m. (37) 

Thus we have established the following result. 
Theorem 2. Assume that curves ^{x,t) on the sphere S"(l) (n > 1) are governed by the flow 

where {ei, 62, • • • , e^} is the natural frame of the curve 7(x,t), {ui,U2, ■ ■ ■ ,Un-i) is defined by 
the curvatures k = fh = u + Uxx, Q = \u\'^ + Then the flow (38) is intrinsic and the 

curvature vector m fulfills the equation (37). 

Remark 2. In the case of n = 2, i.e., the case of S^(l), equation (37) reduces to the scaler 
modified CH equation (4) with 5 = —1, which is completely integrable. In the case of n = 3, let 
ui = u, U2 = V, m = u + Uxx, n = V + Vxx, then system (37) reduces to 

mt = ^[{u^ + v'^ + ul + vl)m]^ - {uvx - vux)n, 

nt = ^ [{u^ + v'^ + ul + v'^)n]^ + {uVx - vux)m. 

In general, the multi-component system (37) can be written as 

^ n— 1 n— 1 

^ 2 ^ f^^^ ^ ■"ix)"^i]^-X]("*^-?>~'"-j'^''^^"^-?'' 1 < « < ■ra-l, mi = Ui+ Ui^xx- 

Case 2. h = Ux,k = rn = u — Uxx- In this case, the tangent velocity / is given by 

1 """^ 1 
/ = (^x^k^ - Uxx,Ux) = 2 (^»^ ~ '^l^) + ci = 2^ + (^^^ 

i=l 

where ci is an integration constant. Substituting (39) with ci = 1 into (35) and noting that 

iyd~^{kx A Kj) k = — {u A Ux) rh — [ux Am)m = — [u A Ux) m — {fh, rh)ux + -^Qxfh, 
we obtain the multi-component modified CH equation 

'mt = ]^{mQ)^- {u Aux)fh + 2ux, m = u-Uxx- (40) 
Thus we have proved the following result. 
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Theorem 3. Assume that curves ^{x,t) on the sphere S"(l) (n > 1) are governed by the flow 

,^ . n-l 

^ / j=l 

where {ei, 62, ... , Cn} is the natural frame of the curves 7(x, t), {ui,U2, ■ ■ ■ , Un-i) is defined by 
the curvatures k = fh = u — Uxx, Q = \u\'^ — \ux\'^- Then the flow (41) is intrinsic and the 
curvature vector fh satisfies (40) . 

Case 3. h = Ux, k = In this case, the tangent velocity / is determined by 

/ = d~^{-ilxx,Ux) = - '^\ux\'^ + C2, 
where C2 is an integration constant. It is inferred from (30) and (31) that 

i= (^{Uxl"^ + 'i^ - C2^ Uxx, 6x = 0. 

Setting = 0, C2 = 1 and v = Ux, we arrive at the multi-component short-pulse equation 

Vxt + l{\v\^Vx)^ + v = 0. (42) 

Hence we could prove the following result. 
Theorem 4. Assume that curves ^{x,t) on the sphere S"(l) (n > 1) are governed by the flow 

n-l 

7i = ( 1 - 7,\ux\^ ] ei + Uj^xej+i, (43) 



where {ei, 62, ... , 6n} is the natural frame of the curves 7(x, t), {ui,U2, ■ ■ ■ , Un-i) is defined by 
the curvatures k = m = —Uxx- Then the flow (43) is intrinsic and the curvature vector v = Ux 
satisfies (42). 

The scalar equation of (42) was derived by Schafer and Wayne [52] as a model for the 
propagation of ultra-short light pulses in silica optical fibers, which is also an approximation of 
nonlinear wave packets in dispersive media in the limit of few cycles on the ultra-short pulse 
scale. 

It is worth noting that, analogous to the system (37), there is another version of multi- 
component generalization of the modified Camassa-Holm equation (4) with 5 = 1 

fht = ^{fhQ)x - {u Aux)'rn, rh = u-Uxx, (44) 

(where Q = |np — juxP). However, in contrast to the system (37) arising from the compact 
Riemannian symmetric space §"(1) = SO(n + 1)/S0(n) with positive constant curvature, the 
system (44) arises from the noncompact Riemannian symmetric space with negative constant 
curvature, i.e., the hyperbolic space H" = SO(n, 1)/S0(n). The derivation to integrable curve 
flows is similar to that for the S'^(l) case. In fact, the Lie algebra structure corresponding to 
the hyperbolic space H" = SO(n, 1)/S0(n) is described by 

so(n, 1) = e m = so(n) © M", 
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with 









j G m, 


\0 





where p G M", G so(n). The Cartan connection matrices (2j{Dx) and u}{Dt) of the natural 
frame are replaced with 



and 





(o 


1 


0^ ^ 




L0{Dx) = 


1 







1 , O G so(n - 1 






A? 








^0 


/ 


/in 




u^iDt) = 









, *Gso(n-l 




yh 









for curves on H" = SO(n, 1)/S0(n). Similar results can be derived for the Hyperbolic space H'^. 
Integrability of the equations (37) and (44) are guaranteed by the following results. 

Theorem 5. The systems (44) and (37) are Lax integrable, namely they admit the following 
(n + 1) X (n + 1) Lax-pair 



y = ( iQ + + ^XQrn^ 

—X^^u—hXQfh —uAUx 



(px = U(f), 




= V(P, 


where for (44), 








1 


^ \ 







Xrn^ 1 




—Am 


/ 


and for (37), 








1 









Xifi^ 




—Xrh 


j 




\Q - A-2 X-^ul 

V=\-\Q + X-^ -X-^u^ + \XQm^' 

—X^^Ux X^^u — \XQrh —uAux 



6 Concluding remarks 



In this paper, geometrical formulations to several multi-component integrable systems are pro- 
vided. These systems are regarded as multi-component generalizations of the CH equation and 
the modified CH equation, which can be obtained through the tri-Hamiltonian duality approach. 
We showed that an integrable generalization to the nonlinear Schrodinger equation arises from 
a non-stretching invariant curve flow in the three-dimensional Euclidean geometry. The in- 
tegrable complex CH equation comes from an invariant curve flow on the Mobius 2-sphere. 
Furthermore, we verified that multi-component generalizations to the modified CH equation 
arise naturally from the curve flows in n-dimensional sphere §""(1) and the hyperbolic space 
H" = SO(n, 1)/S0(n). 

In [27], Olver, Kamran and Tenenblat have established the theory for curves in affine sym- 
plectic geometry. The curve flows in four-dimensional affine symplectic geometry were studied 
in [58] , and an integrable three-component equation with bi-Hamiltonian structure was obtained. 
The theory for curves in the centro-equiaffine symplectic geometry was established in [54]. It 
was shown that certain invariant curve flows in the centro-equiaffine symplectic geometry yield 
noncommutative KdV equations [47] . It is still not clear that what are the dual version of these 
integrable equations arising from curve flows in the affine and centro-equiaffine symplectic 
geometries. 
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